Radio-Frequency Single-Electron Refrigerator 
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We propose a cyclic refrigeration principle based on mesoscopic electron transport. Synchronous 
sequential tunnelling of electrons in a Coulomb-blockaded device, a normal metal-superconductor 
single-electron box, results in a cooling power of ~ &bT x / at temperature T over a wide range 
of cycle frequencies /. Electrostatic work, done by the gate voltage source, removes heat from 
the Coulomb island with an efficiency of ~ ksT/A, where A is the superconducting gap. The 
performance is not affected significantly by non-idealities, for instance by offset charges. We propose 
ways of characterizing the system and of its practical implementation. 



Cyclic conversion of heat into mechanical work, or 
conversely, work into extracted heat, demonstrates basic 
concepts of thermodynamics. The latter forms the basis 
of refrigeration (see, e.g., pj). Whereas cryogenic refrig- 
eration dates back to the 19th century 0, U, electronic 
microrefrigeration has existed only for a little longer than 
a decade (3j- Successful solid-state refrigeration schemes 
operating at cryogenic temperatures, typically in the sub- 
kelvin regime, rely most commonly on normal metal- 
insulator-superconductor (NIS) tunnel junctions 0, Q. 
They have allowed substantial electron 0, as well as 
phonon 0, @ temperature reduction. Quasiparticle re- 
frigeration occurs thanks to the superconducting energy 
gap which allows only most energetic electrons to escape 
in a tunneling process, thus effectively cooling the N re- 
gion. Until now no electronic cyclic refrigerators have, 
to the best of our knowledge, been demonstrated. Yet 
there exists a proposal to use a mesoscopic semiconduc- 
tor ratchet as a Brownian heat engine for electrons 0- 
In this Letter we show that a simple mesoscopic hybrid 
nano-structure in the Coulomb blockade regime exhibits 
significant cooling power when work is done on it by a 
radio-frequency gate, even in the absence of DC-bias. 
The presented concept bears resemblance to that of a 
single-electron pump, where an electron is transported 
synchronously at an operation frequency / through a de- 
vice producing average electric current of magnitude ef 
[lot lilt fl2| . In the present system, heat at temperature 
T is transported similarly due to cyclic gate operation in 
the Coulomb blockade regime with average "current" of 
order k^Tf. 

As a basic but representative example we study a 
single-electron box (SEB) with a normal metal island and 
a superconducting lead [see Fig. H( a )L a system whose 
electrical properties are determined by single-electron 
charging effects p|. In practise it can be easily real- 
ized with present-day fabrication technologies. In SEB, 
a tunnel junction of capacitance C and resistance i?x 
connects the island to the lead. A source providing volt- 
age V g is furthermore capacitively coupled to the island 
through a capacitance C g . The charging energy of the de- 
vice can be written as £7 c h(n, rig) = Eq{u + n g ) 2 . Here, 
Eq = e 2 /(2Cs) is the unit of charging energy, n is the 
number of extra electrons on the island, and n g = C g V g /e 



(a) r ^ r + < b > 



/? T , C_L 






I 


s 


L 






-2 -1 


) 1 2 


r 








g 





(d) 5 

UJ° 3 
■~~c 2 
UJ° 1 






r ■ 


• 













|2 
1 

|o c a 

-1 
-2 



0.0 0.2 0.4 0.6 0.8 1.0 

oi/2ji 



FIG. 1: (Color online) Single-electron refrigerator (SER). (a) 
Single-electron box with a normal metal (N) island and a su- 
perconducting (S) lead, (b) The trajectory on the (n, n g ) 
plane for a representative cycle studied, (c) Sketch of energy 
band diagrams of the SER showing the tunneling processes 
in this cycle, (d) The charging energy of the system (solid 
line, left scale), where discontinuities are observed as elec- 
trons tunnel. This change of energy, provided by the voltage 
source, is supplied to the tunneling electrons to overcome the 
superconducting gap. The gate cycle is shown by the dashed 
line (right scale). 



is the amount of charge in units of e induced by the 
gate on the island. The total capacitance of the island is 
Cs = C + C g + Cq, where Co is the self-capacitance of 
the N region. 

First we study charge and energy transport in the de- 
vice, which allows us to obtain the main result of this 
paper, i.e., heat extraction in a cyclic operation of the 
gate voltage. In the following this system will be referred 
to as single electron refrigerator (SER). We identify two 
stochastic events, tunneling into the island from the lead, 
which we associate with a "+" index, and tunneling out 
from the island to the lead ("— " index). The respective 
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changes in charging energy read: 

e+ = 6E+JA = 2e c {n + n g + 1/2) 

e- = SE-J A = -2e c (n + n g - 1/2). (1) 

Here, A is the superconducting energy gap and ec = 
E c /A. 

Within the theory of sequential single-electron tunnel- 
ing ^3], assuming that the tunneling electron does not 
exchange energy with the environment [13( , the two rates 
of electron transport can be written as 

r+ ( £+ ) = ^ / defl sW/s( e )[l - Me - 6+)] 

r ~(0 = ^ / de " s ( e ^ N ( £ + e ~)[l - /s(e)].(2) 

Here, the normalized BCS density of states at energy E = 
eA reads ns(e) — e/y e 2 — 1. To take into account non- 
idealities in realistic superconductors [l5j , we use ns(e) = 
|Re[(e + ry)/y/{e + i'j) 2 — 1]| instead, where the typical 
value of the smearing parameter 7 is in the range 1 • 
1CP 5 ...1TCP 3 for aluminum as a thin-film superconductor 
0. Furthermore, /s(n)( £ ) is the Fermi-Dirac distribution 
function in S (N) at temperature Ts(T) and energy E. 
The corresponding heat fluxes associated with the two 
tunneling processes are given by 

= ^ / *( c ~ e + )n s (e)/ s (e)[l - / N (e - e+)] 

Q~(n = ^j de(e + e-)n s (e)/ N (e + e-)[l-/s(e)],(3) 

where the first one is the energy deposition rate by incom- 
ing electrons, and the second one is the energy extraction 
rate by leaving electrons. 

We start by an approximate calculation which yields 
illustrative yet nearly quantitative results. At low tem- 
peratures, T, Ts <C A/fce, and in the domain 1 + > 0, 
Eqs. © yield 

J*&*^J^e-&^\ (4) 
With the same approximations, we get from Eqs. J2J 

(Q ± ( e ± ))c T [^ + A(l + e ± )] (5) 

for the average energy per event of a tunneling elec- 
tron with respect to the Fermi energy, (Q ± (e ± )) = 
Q ± (e ± )/r ± (e ± ). We expect the exponent in Eq. (@J 
to assume values in the range 1...10 at that point in the 
cycle where tunneling is most likely to occur. This is be- 
cause the prefactor in front of the exponential term, the 
" attempt frequency" , is typically a few orders of magni- 
tude higher than the operation frequency of the device 
[l6j . The tunneling electron is then expected to have 
energy 

Q ± ~ Tk B T. (6) 




FIG. 2: (Color online) Heat extracted from the N island in 
a SER during the cyclic operation of the gate. These calcu- 
lations were performed with the following parameter values: 
k B T/A = 0.05, e c = E c /A = 0.3, R T = 30 kfi, 7 = 1 ■ 10~ 4 , 
/ = u/(2ir) = 10 MHz, and A = 200 ^eV, which corresponds 
to aluminum as the superconductor. The gray line shows the 
the stochastic results, and the rising black line is the solution 
of the master equation. Sinusoidal line indicates the instan- 
taneous gate position n g , and the piece- wise constant red line 
shows the number of extra electrons on the island in a typical 
trajectory. The inset displays the average number of charges, 
n avc , at each instant within the cycle obtained from the mas- 
ter equation. 

This states that the "+" process brings a negative 
amount of energy into the island, whereas the "— " pro- 
cess removes a positive energy from the island, i.e., all 
the events tend to extract hot excitations from the N is- 
land, which results in refrigeration of it. The expectation 
value of energy carried out from the island in a full cycle 
is then Q — n t (Q~ — Q + ), where nt is the number of tun- 
neling events in both directions within a full cycle. We 
then have Q = rk^T, where r is of the order of n t . By 
repeating cycles at frequency / one obtains an average 
cooling power 

Q = rk B Tf. (7) 

To go beyond these approximate results, we employ ei- 
ther deterministic, master equation-based computations, 
or stochastic simulations of charge and energy dynamics. 
In the first method, we follow the probabilities p n (t) of 
observing charge number n on the island at time instant t. 
Equations J5J) yield for each n state the time-dependent 
transition rates corresponding to the "+" and "— " 
processes. The master equation then reads 

Pn = -(r+ + T~)p n + r+_ x p„_i + r~ +1 p n+1 . (8) 

Furthermore, we write the expectation value of instanta- 
neous heat current with the help of Eqs. © as 

P = Y J Pn{Qn~Qt)- (9) 

n 

We keep track of only a small number (~ k^T/Ec) of 
lowest-energy charge number states, and assume the rest 
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to have zero occupation probability. With a reasonable 
initial condition, e.g., p n (0) — S n o, p n (t) converges to 
an approximate periodic equilibrium during the first few 
cycles. The total cooling power Q is then obtained by 
averaging Eq. © over one cycle. 

In the stochastic method, we employ a Monte Carlo 
simulation of the tunnel events to generate a sample tra- 
jectory in n space. We evaluate the integrands in Eqs. 
(J2J over a small energy interval 8e. Multiplied by a small 
time step St, these yield the probability of observing a 
"+" event at energy e — e + or a "— " event at energy 
e + e~ during St. If the "— " ("+") process occurs, the 
energy is added to (subtracted from) the cooling budget. 
The total cooling power can be obtained by averaging 
over a trajectory long enough to suppress the stochastic 
variation to a negligible level. 

As an example we consider a cycle where the gate volt- 
age varies periodically such that n g = 2sin(a;t), see Fig. 
^ (d). As presented in Fig. [T] (b), n then tends to fol- 
low the sequence +1— » — ► — 1— > — ► +1. The thick 
vertical lines indicate the stochastic uncertainty of the 
tunnel events along the gate excursion. In Fig. ^ (c) 
and (d) we show, respectively, the tunnel events on an 
energy diagram during the example cycle demonstrating 
the cooling feature of the device, and the corresponding 
charging energy along the cycle. Figure shows numer- 
ical results over 15 cycles, of charges on the island and 
heat extracted from it. In the stochastic simulation the 
ideal cycle of n : +1 ^ —1^0-^ +1 is repeated 
almost regularly over these cycles, although the system 
visits once the n — — 2 state and the exact time instants 
of tunneling jitter around their mean values. The deter- 
ministic calculation yields also the expectation value of 
charge on the island at each value of gate (inset). The 
efficiency of the refrigerator, defined as the ratio of heat 
extracted and work done by the source, is r\ ~ 4.0/cbT/A 
in this particular case. This exceeds by more than a 
factor of three the corresponding figure of a static NIS 
refrigerator, rj ~ 1.2/cbT/A l 3]. This difference stems 
from yet another favourable energy filtering in SER: due 
to Coulomb blockade, only one electron tunnels in each 
phase of the cycle, and it extracts maximal amount of 
heat from N. 

Further insight into the performance of SER can be 
gained by inspecting the frequency dependence of Q, 
shown in Fig. (a) for different values of Rt- Q is a 
non-monotonic function of /, decreasing both at low and 
high values of /. Increasing the junction transparency 
leads to an enhancement of the maximum heat current 
and to widening of the window of operational frequen- 
cies. This can be easily understood in terms of reduced 
discharging time through the junction. For Rt = 30 kSl 
(5 kQ), Q is maximized around / ~ 150 MHz (1 GHz), 
where it obtains values of ~ 1 fW (~ 6 fW) for Al. As 
will be shown below, such values will lead to substantial 
electronic temperature reduction in the N island. 

Whether a single-electron device is usable in prac- 
tice depends on its immunity to fluctuating background 
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FIG. 3: (Color online) (a) Heat flux extracted from the N is- 
land of a SER as a function of frequency for some values of Rt. 
From bottom to top: Rt = 100, 30, 10, and 5 kfi. The compu- 
tations were performed assuming n s (t) = 2sin(u;£). (b) Nor- 
malized energy extracted per cycle versus background charge 
no. In the calculations we used / = 10 MHz, and Rt = 30 
kSl. The gate amplitudes for the three curves from bottom 
to top are 2, 3.5, and 5. In both (a) and (b) we assumed 
A = 200 A*eV, k B T/A = 0.05, E c /A = 0.3, and 7 = 1 ■ 10~ 4 . 



charges (Qa). Their influence can be studied numerically 
by adding an offset no = Qo/e to n s . The results are dis- 
played in Fig. 13(b), which indicates heat extracted per 
cycle vs no for three values of gate amplitude. Clearly no 
has some effect on the cooling power of SER, but these 
small variations can be compensated for by adjusting ei- 
ther the amplitude or the DC-bias of the gate, if neces- 
sary. Furthermore, in practice it may be advantageous 
to suppress the influence of no by employing larger gate 
amplitudes at lower /, which leads to the same cooling 
power, but to weaker no dependence. 

In order to test the performance of the SER in a realis- 
tic context, we need to take into account the mechanisms 
that drive power into the N island. We expect the sys- 
tem to be in the quasi-equilibrium limit |3j , where strong 
inelastic electron-electron interaction forces the quasipar- 
ticle distribution toward a thermal one at a temperature 
T which can differ from that of the lattice (Tbath)- At 
low temperatures, which for metallic structures are typ- 
ically below 1 K, the main heat load on electrons comes 
from phonon scattering with energy exchange according 
to Q c -bath(T,T bath ) = £V(T 5 - T b 5 ath ) 0, where E is 
the electron-phonon coupling constant and V is the is- 
land volume. In order to cool the island, frequency / 
has to be high enough to prevent full relaxation toward 
Tbath- Thus we require (r _p h ) _1 -C / < (t^)^, where 
r e-ph(e-e) is the electron-phonon (electron-electron) inelas- 
tic scattering time. These constraints can be met exper- 
imentally recalling that, e.g., for a copper (Cu) island 
one has (r c _ph) _1 < 10 5 Hz at T < 200 mK 18], while 
(r e _e) _1 ~ 10 9 — 10 10 Hz even in quite pure metal sam- 
ples |l9j . Heat input by other relaxation mechanisms, 
in particular radiative electron-photon heat load, can be 
ignored in this case because of poor matching [2£J, UM ■ 

In Fig. 4 (a) we plot the steady state temperature of an 
N island of a typical size. The results are shown as a func- 
tion of / at few representative values of Tb a th, which we 
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FIG. 4: (Color online) (a) Temperature of the N island vs 
gate frequency at a few bath temperatures. We assume that 
the volume of the island is V = 1 • 10" 21 m 3 , E = 1 • 10 9 
WK- 5 m" 3 , £ C /A = 0.3, and 7 = MO -4 . T c = A/(1.764fc B ) 
is the superconducting critical temperature. The results are 
shown for Rt = 10 ktt (solid lines) and Rt = 30 ktt (dashed 
lines), and for bath temperatures (from bottom to top) 131 
mK (T bath /T c = 0.1, black), 197 mK (T bath /T c = 0.15, red), 
263 mK (r ba t h /rc = 0.2, blue), 328 mK (T bath /T c = 0.25, 
green) and 394 mK (Tb at h/Tc = 0.3, brown), (b) A scheme 
for refrigerating a device on a thin dielectric platform using an 
array of SERs. The SERs cool down electrons and phonons 
on the membrane, as well as the device on the platform. The 
islands share a common gate and a superconducting lead out- 
side the membrane. 



furthermore assume to be equal to Tg. The lines are solu- 



tions of the balance equation Q(T) + Q e -b&th{T, Tb a th) = 
0, where Q(T) is the cooling power of SER at tempera- 
ture T. The results demonstrate significant temperature 
reduction over a broad range of frequencies. In practise 
T can be measured, e.g., by an additional NIS junction, 
which probes the temperature dependence of a tiny quasi- 
particle current. A thermometer based on supercurrent 
in a Coulomb-blockaded SINIS junction is another possi- 
bility 22] . In Fig. 4 (b) we show a possible implementa- 
tion of an array of SERs to refrigerate a separate device, 
e.g., a radiation detector on a dielectric membrane In 
this realization, the performance of SERs can be probed 
by measuring the drop of lattice temperature ■ 

We presented the concept and analyzed the perfor- 
mance of a cyclic single-electron refrigerator, which in- 
troduces thermodynamic cycles in a nano-electronic sys- 
tem. The concept can readily be generalized to more 
complex structures with further enhanced performance. 
For example, a three terminal device, in form of a SINIS 
single-electron transistor, can be controlled additionally 
by applying a drain-to-source voltage across. Coulomb 
blockade suppresses the cooling power of a statically bi- 
ased SINIS refrigerator [23(, but adding radio- frequency 
gating can enhance the cooling power even above that 
of the SER presented here. Finally, we analyzed here 
quantitatively only sinusoidal cycles of n g in the basic 
SER. Alternative architectures and cycles may lead to 
better control between Q and ksTf, and the device might 
then yield a "quantized" and synchronized heat current 
in analogy to electrical current in metrological devices. 
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